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ABSTRACT 



In this fourth part of the series presenting the Optical Gravitational Lensing Ex- 
periment (OGLE) microlensing studies of the dark matter halo compact objects (MA- 
CHOs) we describe results of the OGLE-III monitoring of the Small Magellanic Cloud 
(SMC). Three sound candidates for microlensing events were found and yielded the 
optical depth tsmc-oiii = 1-30 ± 1.01 x 10~ 7 , consistent with the expected contri- 
bution from Galactic disk and SMC self-lensing. We report that event OGLE-SMC-03 
is the most likely a thick disk lens candidate, the first of such type found towards the 
SMC. In this paper we also combined all OGLE Large and Small Magellanic Cloud 
microlensing results in order to refine the conclusions on MACHOs. All but one of 
OGLE events are most likely caused by the lensing by known populations of stars, 
therefore we concluded that there is no need for introducing any special dark matter 
compact objects in order to explain the observed events rates. Potential black hole 
event indicates that similar lenses can contribute only about 2 per cent to the total 
mass of the halo, which is still in agreement with the expected number of such objects. 

Key words: Cosmology: Dark Matter, Gravitational Lensing, Galaxy: Structure, 
Halo, Galaxies: Small Magellanic Cloud 



1 INTRODUCTION 

Campaigns like OGLE, MACHO or EROS were all founded 
nearly 20 years ago primarily to settle the question whether 
dark matter resided in the galactic halo in the form of 
compact objects (MACHOs). The background theory of the 
method employing gravitational microlensing was provided 
by Bohdan Paczyhski (Paczyhski 1986) and in subsequent 
years Paczyhski himself advocated and supported the cre- 
ation of large-scale photometric surveys towards the Magel- 
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lanic Clouds. The Large and Small Magellanic Clouds (LMC 
and SMC) were considered the best targets for such studies, 
as with a moderate telescope size the stars of the Clouds 
were easily resolved and hence formed a perfect background 
for the hunt for microlensing events. These are unique and 
temporal brightening events observed on a source back- 
ground star caused by the gravitational lensing phenomenon 
in which the lensing object is located in front of the back- 
ground star. 

In the most generic case of a point lens mass and a 
point source the gravitational amplification of the back- 
ground light is described by equation (after Paczyhski 1996): 
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where to is the time of the maximum of the peak, £e is the 
Einstein radius crossing time (event time-scale) and uo is 
the event impact parameter. 

Parameter £e is the only parameter linked with the 
physical properties of the source-lens system and is depen- 
dent on the lens mass, the lens and source distances and the 
relative velocity between lens and source. This means that, 
in principle, a single microlensing event can not provide any 
information regarding the lensing object. Additional data 
are possible only in case of non-standard events, like bi- 
nary lens/source or those showing a parallax effect. Hence, in 
studying the mass distribution of lenses, microlensing events 
are regarded only in groups. The statistical value related to 
the total mass enclosed in the volume towards the back- 
ground sources is the microlensing optical depth, r: 

Afov , 

r = - — (2) 

2A/»T obs 4- e(tEi) { ' 

where T b s is the time-span of all observations, N* is the 
total number of monitored stars, N cv is the total number 
of events, each with time-scale £e, and detection efficiency 
e(*Ei). 

In the case of the Magellanic Clouds the overall opti- 
cal depth can be composed with contributions from lenses 
at various levels. Apart from the hypothetical dark matter 
compact objects there can also be a contribution to the opti- 
cal depth from luminous lenses from the disk of our Galaxy 
and from the frontal side of the Cloud itself (i.e., self- lensing, 
SL). 

Observations and studies of the pixel-lensing microlens- 
ing events towards M31 are still not conclusive with respect 
to the amount of self-lensing and MACHO lenses (see Calchi 
Novati 2010b for recent review). However, the most recent 
analysis of the data tend to be inclined towards self-lensing 
as being the only explanation of the observed signal (Calchi 
Novati 2010a). 

The most recent studies of the data gathered by mi- 
crolensing surveys towards the LMC suggest the contribu- 
tion from dark matter compact objects with masses below 
10 Mq can be practically ruled out. The upper limit on 
MACHO abundance in the halo was set to 6-7 per cent 
by the EROS (Tisserand et al. 2007) and OGLE groups 
(Wyrzykowski et al. 2011) (see Moniez 2010 for review). 
The core of the detected microlensing signal seems to come 
from the self-lensing events, as suggested already in Sahu 
(1994b) and Sahu (1994a). 

The LMC self-lensing effect is relatively well understood 
thanks to the convenient, almost face-on alignment of the 
LMC towards us. The studies by Mancini et al. (2004), 
Calchi Novati et al. (2009) and Calchi Novati & Mancini 
(2011) showed that the average self-lensing optical depth is 
in the order of 0.4 x 10~ 7 , which is consistent with what 
was measured in OGLE-II (Wyrzykowski et al. 2009, here- 
after Paper I) and OGLE-III data (Wyrzykowski et al. 2011, 
hereafter Paper III). 

The situation with the SMC is more complicated. 
This dwarf galaxy seems to be elongated and is stretched 
out along the line-of-sight, which produces a much higher 
self-lensing optical depth. Some estimates by Palanque- 
Delabrouille et al. (1998) and Graff & Gardiner (1999) place 
this value at at least 1.0 x 10~ 7 , but it can actually be much 



higher. Microlensing studies of the SMC tend to confirm this 
with their typically higher event rates (Tisserand et al. 2007, 
Wyrzykowski et al. 2010, hereafter Paper II), compared with 
the LMC. 

In this paper we conclude the series of studies of 
the OGLE microlensing data towards the both Magellanic 
Clouds with the analysis of the 8 years of observations of the 
SMC in the course of the OGLE-III project. We attempt to 
finalise the investigation of the subtle topic of dark matter 
compact halo objects, providing the result from currently 
the best and the most suitable data set. 

The paper is organised as follows. First, the OGLE- 
III data used for the analysis are presented. Then, we de- 
scribe the method applied for finding the microlensing events 
and its results with a detailed description and study of each 
found event. This is followed by a calculation of the optical 
depth and discussion of the results. 



2 OBSERVATIONAL DATA 

The photometric data used in this study were collected dur- 
ing the third phase of the OGLE project (2001-2009) with 
the 1.3-m Warsaw telescope located at Las Campanas Ob- 
servatory, Chile, operated by the Carnegie Institution of 
Washington. The "second generation" camera comprised 
eight SITe 2048 x 4096 CCD detectors with 15 /im pixels 
resulted in 0.26 arcsec/pixel and 35 x 35 arcmins total field 
of view. The details on the instrumentation setup can be 
found in (Udalski 2003). 

The central regions and the outskirts of the SMC were 
covered by 41 fields, giving a total of 14 square degrees. The 
map of the fields is shown in Fig. 1. The statistical details of 
the fields are gathered in Table 1 with the following informa- 
tion: the coordinates of their centers, the number of "good" 
template objects in the /-band, the blending-corrected num- 
ber of stars (see Section 5) and the mean number of all 
objects visible on a single CCD (of 8) used for assigning 
the blending density level. "Good" objects are the template 
objects with at least 80 good data points (excluding mea- 
surements with very large error-bars) and mean magnitude 
brighter than 21.0 mag (chosen as a mean peak of the ob- 
served luminosity functions). 

The very first observations of the SMC within the 
OGLE-III phase were taken in June 2001 (JD=2 452 085), 
except field SMC 140 (started in July 2004), and continued 
until May 2009 (JD=2 454 954). Most of the observations 
were taken through the Cousins /-band filter with exposure 
time 180 s with a mean seeing of 1.36 arcsec. Between 583 
and 762 measurements were gathered in each field with an 
average sampling varying from 2.0 to 3.0 days between sub- 
sequent frames (excluding the gaps between the seasons). 
The only exception is field SMC128 which was observed 1228 
times with an average sampling of 1.6 days. Additionally, 
between 47 and 114 observations per field were obtained in 
Johnson V-band with integration time 225 s and a mean 
seeing of 1.39 arcsec. First V-band observations were taken 
in November 2004 (JD=2 453 314) and since then V filter 
frames were taken regularly. The average sampling frequency 
in the V-band was between 4.6 and 10.3 days. 

The image reduction pipeline used an image subtraction 
technique and was based on Difference Image Analysis (DIA; 
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Table 2. Error correction coefficients for each CCD chip of the 
first four OGLE-III LMC fields for I- and V— bands. The full 
table is available on-line from the OGLE website. 
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.145 


0.0044 


1 
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0.0026 
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.951 
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Alard & Lupton 1998, Wozniak 2000). The data used in this 
work come from the final reductions calibrated to the stan- 
dard system. The template images used for subtraction were 
composed from selected best individual frames resulting in 
a seeing of around 0.9 arcsec and 1.0 arcsec for I-band and 
V-band, respectively. The full description of the reduction 
techniques, photometric calibration and astrometric trans- 
formations can be found in (Udalski et al. 2008). 

Photometric errors produced by the difference imaging 
process were adjusted as described in Paper I by deriving 
a magnitude dependent correction factor using the constant 
stars. For each SMC field and each filter (/ = I, V) we found 
parameters 7/ and e/, which correct the original error-bar 
according to the formula: 

Omagf cor = \J (7/°'niagf ) 2 + e p (3) 

where (J mag is the original error-bar returned by the DIA. 

Error-correction parameters averaged over all SMC 
fields yielded with (7/) = 1.04031, (ej) = 0.00414852, 
(7v> = 0.867849, (t v ) = 0.0026929. Table 2 shows the error 
correction coefficients for the first few fields. The full table 
can be found on the OGLE website 1 . 



3 SEARCH PROCEDURE 

For the SMC OGLE-III data we applied the same auto- 
mated pipeline for detecting microlensing events as designed 
and used in searching through the LMC OGLE-II (Paper I), 
SMC OGLE-II (Paper II) and LMC OGLE-III (paper III) 
data. The procedure for the SMC contains 10 cuts applied 
to the pre-selected "good" data (with enough good obser- 
vations) and limited to the right magnitude level. Table 3 
describes all cuts applied to the data for two star samples 

1 http://ogle.astrouw.edu.pl/ 



(defined at cut 0): All Stars - limited to 21 mag (maximum 
of the observed luminosity functions, i.e., the end of OGLE 
completeness in these fields) and Bright Stars - limited to 
19.3 mag (1 mag below mean Red Clump position). 

In the first cut for each star sample we selected light 
curves with a significant bump over baseline, following the 
concept of Sumi et al. (2006). Then, we removed all candi- 
dates with their baseline colour and magnitude located in 
the "blue bumper" region to exclude these most common 
contaminants related with the Be stars. In cut 3 we chose 
those light curves for which the standard microlensing fit 
(with blending parameter fixed, hence index /Lt4 in the table 
3) was better than a constant line fit. In the following cut 
we required that there were at least 6 data points in the 
microlensing peak between to — te and to + £e- 

In cut 5 we fitted the candidate bump also with the 
simple supernova model (exponential decline) and selected 
only those for which the \ 2 of the microlensing fit was better 
(either with blending or without). The number of expected 
supernovae in the entire OGLE-III coverage of the SMC is 
around 22 (after Dilday et al. 2010, assuming 20 per cent de- 
tection efficiency) , and indeed many of the objects removed 
at this stage resembled supernovae outbursts, some clearly 
with the host galaxy visible on the OGLE image. In this way 
we also removed all other asymmetric bumps. The objects 
remaining after this step were visually inspected and most of 
them turned out to be either nova-like events or red-baseline 
bumpers. 

We then checked if the peak of the candidate bump 
was located within the data span (cut 6). In the next two 
cuts (7 and 8) we required that the blended microlensing 
model fit converged and yielded reasonable values, allowing 
for some negative blending (effect of instrumental origin, see 
e.g., Smith et al. 2007), and that the \ 2 /N do f of the blended 
fit was better than 2.6 and for non-blended fit around the 
peak better than 4.5. 

Finally, we restricted our search of microlensing events 
to the ones with time-scales between 1 and 1000 days and 
with impact parameters below 1. 



4 RESULTS 

Running our search pipeline on the SMC OGLE-III data 
found 3 candidates for microlensing events in All Stars Sam- 
ple. Two of them were found in the Bright Stars Sample. The 
three events were dubbed OGLE-SMC-02, OGLE-SMC-03 
and OGLE-SMC-04, with the numbers following the num- 
bering started with the single candidate detected in the 
OGLE-II SMC data (Paper II). 

The basic information about these candidates are shown 
in Table 4. First two events were detected while on-going by 
the Early Warning System (EWS) of OGLE, while the third 
one was not known before and occurred at the beginning of 
the OGLE-III when EWS was not yet operating. 

Figure 2 shows the colour-magnitude diagram with can- 
didate events found in OGLE-II, OGLE-III and MACHO 
data. Figure 3 presents a density of Red Clump stars of 
the SMC with marked positions of all candidate events. 
The light curves of three OGLE-III events with microlens- 
ing models fits are shown in Figure 4, whereas their finding 
charts are presented in Figure 5. 



4 L. Wyrzykowski et al. 




Figure 1. Positions of the OGLE-III SMC fields (black). Also shown are all OGLE-fl fields (red rectangles). The three small filled 
squares show the positions of the HST fields used for our blending determination. Background image credit: ASAS all sky survey. 

Table 1. OGLE-HI SMC fields. 
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19 


23 


89 


9 


94.0 
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351.3 


4 


68 


SMC105 
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41 
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66 
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17 
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4 


66 
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15 
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total 5470352 5971776 



Note: Coordinates point to the centre of the field (centre of the mosaic), each being 35' X 35'. Number of "good" objects in the 
template is provided (N > 80 and (/} < 21.0 mag) together with the estimated number of real monitored stars (see Section 5). Mean 
number of all objects detected on a single CCD used for calculating the density of a field is given in the last column. 



Below we discuss each candidate in detail. 4.1 OGLE-SMC-02 



This event was already known before this work as it was de- 
tected in real-time with the OGLE's Early Warning System 
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Table 3. Selection criteria for the search for microlcnsing events in the OGLE-III SMC data and the number of objects left after each 
cut for the All Stars and Bright Stars Samples. 



Cut no. 



No. of objects left 
All Bright 



Oa 
Ob 



Selection of "good" objects 
Significant bump over baseline 
"Bumper" cutt 

Microlcnsing fit better than constant line fit 

Number of points at the peak* 

Microlcnsing fit better than supernova fit, 

Peak within the data span 
[HJD-2450000] 

Blended fit converged 



N > 80, (I) sC 21.0 mag 5,470,352 
N > 80, (J) sC 19.3 mag 



<?i > 30.0 

peak 

(/) > 19.0 mag, (V -I) > 0.5 mag 
*Le-x^4 > 100 



Npeak > 5 

X | JV >MIN(x 2 ,x^ 4 ) 
2085 < t ^ 4954 

< fs < 1-4 



1,217 

482 
186 

182 
115 
112 

35 



1,614,327 
1,060 

380 
107 

104 

68 
65 

23 



10 



Conditions on goodness of microlcnsing fit 
(global and at the peak) 

Time-scale cut 
[d] 

Impact parameter cut 



s= 2.6 and ■■ X ' j4 - pcak s= 4.5 



1 ^ f E ^ 1000 
< «o C 1 



t magnitudes as in the field SMC100.1 (shifted according to the position of the center of Red Clump) 
*in the range of to ± He 

Table 4. Microlensing events candidates detected in the OGLE-III SMC data. Events detected by the OGLE's Early Warning System 
(EWS) have their EWS designation given in the name column. 



Event's name 
(EWS) 



RA Dec 
[J2000.0] [J2000.0] 



field db baseline / baseline V 

star id [mag] [mag] 



source / source V — I 



[mag] 



lag] 



OGLE-SMC-02 0:40:28.10 -73:44:46.5 SMC128.5 15192 18.368±0.002 19.389±0.004 18.368±0.002* 1.022±0.005* 
(2005-SMC-001) 



OGLE-SMC-03 1:02:37.50 -73:22:25.9 SMC111.7 18094 20.66±0.01 20.77±0.01 20.43±0.32 
(2008-SMC-001) 



O.lOiO.01 



OGLE-SMC-04 1:00:59.32 -73:14:17.1 SMC106.3 33618 18.477±0.001 19.463±0.005 18.69±0.79 0.75±0.10t 



Source brightness and colour was set to this of the event's baseline, assuming no blending and dark lens (after Dong et al. 2007) 

t Source colour was derived using EROS B-band data. 



(EWS) and was designated 2005-SMC-001 2 . It was proba- 
bly one of the longest, the brightest and the best observed 
events found to date in the entire history of microlensing 



http: / /ogle. astrouw.edu. pi /ogle3/ews/2005 /smc-001 .html 



monitoring of both Magellanic Clouds. Its light curve ex- 
hibited not only a clear parallax effect but also a devia- 
tion to the single microlensing model around its peak. This 
anomaly is so tiny that it is not clearly seen in the OGLE 
data alone as shown in Fig. 4. However, its presence was ex- 
plained with a non-caustic crossing binary lens model by 
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Table 5. Parameters of the standard Paczyriski microlcnsing model fits to the OGLE-III SMC events. 



OGLE-SMC-02 

parameter 5-parametcr fit 4-parameter fit 7-parameter fit 



H) 


ooyo.u 




3593.7 


±0.02 


ooyo.u 


+n 09 


t B 


195.6 


±1.9 


163.2 


±0.3 


190.6 


±1.6 


wo 


0.07132 


±0.08208 


0.08915 ±0.00016 


0.07340 


±0.00071 


Io 


18.37 


±0.002 


18.36 


±0.002 


18.37 


±0.002 


fSi 


0.8024 


±0.0089 


1.0 




0.8248 


±0.0078 


V 










19.39 


±0.004 


/s v 










0.9102 


±0.0095 


x 2 


2836.7 




3259.0 




3593.0 




N do! 


2.32 




2.66 




2.62 








OGLE-SMC-03 








parameter 


5-parameter fit 


4-parameter fit 


7-parameter fit 


tn 

H) 


4476.1 


±0.1 


4476.1 


±0.1 


4476.1 


±0.1 


t B 


45.5 


±6.2 


54.1 


±0.9 


47.1 


±4.4 


no 


0.147 


±0.024 


0.1161 


±0.0013 


0.140 


±0.016 


Io 


20.62 


±0.01 


20.62 


±0.01 


20.63 


±0.01 


fSj 


1.32 


±0.26 


1.0 




1.25 


±0.16 


V 










20.77 


±0.01 


/s v 










1.20 


±0.16 


x 2 


708.49 




710.65 




784.13 




2 

AT dof 


1.01 




1.01 




1.02 








OGLE-SMC-04 








parameter 


5-parametcr fit 


4-parameter fit 


7-paramcter fit* 


to 


2611.6 


±0.1 


2611.6 


±0.1 


2611.6 


±0.1 


*E 


18.60 


+1.96 

-1.85 


17.17 


±0.27 


18.3 


±1.8 


UQ 


0.3143 


+0.0595 
-0.0472 


0.3589 


±0.0039 


0.3231 


±0.050 


Io 


18.477 


±0.001 


18.477 


±0.001 


18.477 


±0.001 


fSi 


0.84 


+0.22 
-0.16 


1.0 




0.87 


+0.21 
-0.16 


B 










19.48 


±0.01 


/s B 










1.03 


+0.25 
-0.18 


x 2 


1052.3 




1052.9 




1154.6 




- 


1.43 




1.43 




1.15 





* fit to the OGLE /-band and EROS B-band data 
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Figure 2. Colour-magnitude diagram showing HST (black dots) and OGLE (green) stars of one of the OGLE-III fields. Red dots show 
estimated colour and brightness of the source in all OGLE events, whereas black open circles indicate the baseline of each event. Blue 
and green dots mark MACHO events from 1997 and 1998. 
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Figure 3. Localisation of all SMC OGLE and MACHO events on the Red Clump stars density map. The colour coding is as that used 
in Figure 2. Contours follow the density of Red Clump stars. 



Dong et al. (2007). In their analysis they also successfully 
employed for the first time space-parallax effect from ad- 
ditional observations taken by the Spitzer satellite, what 
helped break model degeneracies and solve the microlens- 
ing event. Dong et al. (2007) concluded that the most likely 
location of the lens is the Galactic halo and with no light 
detected from the lens they suggested the lens is a binary 



black hole with a total mass of around 10M© . Nevertheless, 
they have not completely excluded the self-lensing scenario 
with both source and the lens residing in the SMC. 

In our study we performed a simple modelling of a single 
point lens event with no parallax included. This was the 
model we used in the automated search pipeline, therefore 
we were able to derive the detection efficiency for this event. 
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Figure 4. Light curves and microlensing models of candidates for microlensing events detected in the OGLE-III SMC data. The standard 
microlensing model best-fitting to the data is shown in solid lines with black points and red curve for /-band, and blue points and green 
curve for V-band data (except OGLE-SMC-04, where it shows EROS B-band data). Dashed line shows parallax model fit. The residuals 
of the model fitting are shown as a respective lines and data points at the bottom of each panel. 




Figure 5. Finding charts of the three candidates for microlensing 
The side of each chart is 26 arc seconds. Charts with events 02, 
which the microlensing brightening was detected. 



events from OGLE-III SMC data. East is to the right, North is down. 
03 and 04 are shown from left to right. A cross marks the object on 



The values of the fitted parameters are gathered in Table 
4 and the fit is shown in Fig. 4. The time-scale of £e = 
190. 6± 1.6 days we obtained for I and V data is conveniently 
in rough agreement with the time-scale derived by Dong 
et al. (2007) (between 160 and 190 days, depending on the 
details of the model configuration). 

The blending parameters for both bands obtained in 
our model are relatively close to 1, however not including 



the parallax and the binarity in our model should affect 
this value severely. Full modelling of this event performed 
by Dong et al. (2007) (their Table 1) returned a blending 
fraction oscillating around zero-blending solutions, even sug- 
gesting some amount of negative blending in case of a few 
models. 

We analysed the astrometry of the residual images of 
this event obtained with the DIA and measured the centroid 
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shift caused by blending during the event with respect to 
its position during the baseline as measured on the superb 
quality template image. The accuracy of such a centroid shift 
can reach a few tens of milli-arc-seconds for bright events. In 
case of OGLE-SMC-02 we did not detect any displacement 
of the centroid of light, which indicates there is no additional 
blending light coming from a coincidental neighbouring star. 
The lack of blending in microlensing models by Dong et al. 
(2007) was also confirmed by them with high resolution HST 
imaging and, along with other factors, led to the conclusion 
the lens in this event is either dark or very faint. 

Cross-match of the lens position with available high en- 
ergy data from ROSAT, XMM (F. Haberl & R. Sturm, pri- 
vate communication) and INTEGRAL (A. Frankowski, pri- 
vate communication) returned no signal which can be asso- 
ciated with the black-hole. On the other hand, we would see 
the X-ray signal only if there was an accreting disk around 
the black-hole, the feature very unlikely to be present in a 
binary black-hole system. Further detailed observations and 
studies are required in order to confirm the black hole origin 
of this event. 



4.2 OGLE-SMC-03 

This event was detected in real-time while it was on-going in 
2008 by the EWS (2008-SMC-001) 3 . It occurred on a faint 
blue star but was quite well covered in both OGLE bands 
(see Fig. 4), especially on the rising part of the light curve 
which assures this was not a nova- like outburst. 

In order to improve the quality of the light curve of 
this faint object the original data were re-reduced so the 
flux variations were measured at the actual position of the 
event. The data were then fit with a standard microlensing 
model and model with parallax effect included as shown in 
Fig. 4. The parameters of the models are gathered in Tables 
5 and 6 for standard and parallax models, respectively. 

All models had a strong indication of negative blending 
most likely caused by an underestimation of the template 
flux of the faint source and fringing present at the template 
image. This is not surprising for the object fainter than 20th 
magnitude where unresolved background stars can cause lo- 
cal overestimation or underestimation of the background 
level. Removal of the lensed source following the method 
of Gould & An (2002) revealed no blended objects visible 
above the background of the images. 

The position of the source of the event shown on the 
CMD (Fig. 2) is based on the colour and brightness of the 
source derived from the best microlensing fit with the neg- 
ative blending, but the uncertainty caused by this effect is 
encoded in the large error-bar of the source brightness. 

Because the event was relatively long (it lasted for about 
100 days) we also investigated if the parallax model can re- 
produce the data better than the standard model. The mod- 
elling with MCMC found four equally valid parallax solution 
covering all possible degeneracies. Their \ 2 f° r multi-band 
data fit, however, differed insignificantly (Ax 2 ~ 0.4, see 
Table 6) and were better by A% 2 ~ 3.4 than the best non- 
parallax model with x 2 — 784.13. 

We adopted the geocentric formalism to describe the 




^ CO CNl 



Figure 6. Likelihoods of the North and East components of the 
parallax signal detected in OGLE-SMC-03 event. Two regions 
corresponding to solutions 1 & 2 and solutions 3 & 4 are created 
through the jerk-parallax microlens degeneracy. The directions 
parallel (|j) and perpendicular (_L) to the Sun's apparent accel- 
eration in the geocentric frame are also shown. (See Figure 3 in 
Gould 2004). Grey-scale contour shows boundaries of a levels, 
which correspond to the ratios of the likelihoods of the given re- 
gion (L) to the most likely region (L max ), as derived with MCMC. 
Black cross mark origin of the coordinate system. 




*eff (= w <e) [days] 



Figure 7. Four microlensing solutions generated by the jerk- 
paralax degeneracy and the tto degeneracy in the plane of t e g 
and 7Te,±- The 7r Ej x is a projection of the parallax vector onto 
direction perpendicular to the Sun's acceleration in the geocentric 



frame at the specified time to,j 



this is the natural coordinate 



for the jerk-parallax degeneracy. We use t e g instead of uq since 
this parameter yields much smaller scatter. Numbers indicate the 
individual solutions from Table 6. 



http: / /ogle. astrouw.edu. pi /ogle3/ews/2005 /smc-001 .html 
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Table 6. Parameters of the microlcnsing models with parallax fit 
to the event OGLE-SMC-03. 



model no 


1 

U() > 
(te,± > 


2 

llQ <C 
-Tj,x/2) 


3 

uq > 

(*E,± < 


4 

Uq < 
-Tj,x/2) 


to 


4476.238 


4476.242 


4476.25 


4476.241 




±0.085 


±0.096 


±0.11 


±0.072 


*E 


47 


54 


56 


41 




±21 


±27 


±28 


±11 


TE,JV 


-0.43 


-0.56 


2.31 


2.55 




±0.38 


±0.40 


±0.37 


±0.39 


^E,E 


0.33 


0.37 


0.35 


0.30 


±0.23 


±0.25 


±0.25 


±0.22 


UO 


0.148 


-0.140 


0.145 


-0.152 




±0.036 


±0.039 


±0.039 


±0.033 


h 


20.62 


20.62 


20.62 


20.62 




±0.01 


±0.01 


±0.01 


±0.01 


fSj 


1.41 


1.35 


1.49 


1.47 


±0.34 


±0.42 


±0.36 


±0.35 


V 


20.77 


20.77 


20.77 


20.77 




±0.01 


±0.01 


±0.01 


±0.01 


fSy 


1.47 


1.39 


1.54 


1.54 


±0.35 


±0.42 


±0.36 


±0.38 


x 2 

x 2 

^dof 


780.80 
1.0180 


780.80 
1.0180 


780.49 
1.0176 


780.39 
1.0175 



parallax effects (An et al. 2002, Gould 2004), where we as- 
sume the reference system to be located at the Earth at 
the time to, pa r = 2454476.23, which we choose close to the 
peak of the event. This ensures that the values of uo, to 
and <e in the parallax model are close to those from non- 
parallax model. Any deviations of the observer position due 
to Earth motion are calculated against this reference system, 
which is in rectilinear motion at the velocity of the Earth at 
that time. The heliocentric velocity of the reference system 
projected on the plane perpendicular to the line of sight is 
approximately 29 kms -1 to South. Usage of the geocentric 
system helps to find and identify all degenerate solutions. 

To find all solutions we investigated the four-fold mi- 
crolensing degeneracy that is created in case of a single lens 
with parallax from the two coexisting discrete degeneracies: 
a two-fold uo degeneracy (Smith et al. 2003) and the jerk- 
parallax degeneracy (Gould 2004). These yield four solutions 
that are presented in Figures 6 and 7 - fit parameters for 
solutions 1 though 4 are gathered in Table 6. 

Pairs of parallax solutions indicate either very close (so- 
lutions 3 & 4) or a bit more distant lens (solutions 1 & 2). 
The former would yield a distance to the lens of order of 
100 pc for a typical mass of the lens ~ 0.3 Mq - it is ex- 
tremely unlikely for the potential lens to be located so close 
to the observer. On the other hand, the latter set of solutions 
have parallax scale of 7te = 0.69 ± 0.32. This yields distance 
2.2 ± 4.3kpc for a typical mass, and since it covers signifi- 



cantly more volume in the Galaxy and many more potential 
lenses we choose these as preferred solutions. The detection 
of the parallax signal in this event implies that the lens most 
probably comes from the Galactic thick disk population of 
stars. 

Assuming the microlensing source to be at the distance 
of the SMC, as suggested by the CMD location of the source, 
we can neglect source motion and calculate lens velocity 
projected on the Earth plane to be AU 1 i^ 1 w 50 kms -1 . 
This, when taking into account the angle of the relative lens 
motion and the velocity of the geocentric reference system, 
leads to heliocentric lens velocity, projected on the plane of 
the observer, to be approximately (—70, 30) kms -1 in North 
and East direction, respectively. Since the lens is much closer 
than the source we can neglect the effect of projection and 
note that the true heliocentric velocity of the lens will be of 
the same order. This type of kinematics is not surprising for 
the thick disk object. 

Because there is no blend visible we can derive an upper 
limit on the lens brightness equal to around 21.5 mag in 
/-band. This limit is met by a late M dwarf (with mass 
< 0.2Mq) at the distance of more than 1.5 kpc, and an 
early M dwarf (with mass < O.3M0) at the distance of more 
than 5 kpc. Although we cannot rule out the lens to be a 
dark object in the Galactic halo the lens as a disk star is 
a much more natural explanation. Calchi Novati & Mancini 
(2011) showed in their simulation of the microlensing rate 
towards the LMC that the contribution of the Galactic Disc 
lenses to the overall microlensing rate is comparable to the 
contribution from the LMC self-lensing. It means that in our 
searches for microlensing events towards Magellanic Clouds 
we should expect some of them to be caused by Galactic 
lenses. 

We can conclude that the most likely scenario for 
OGLE-SMC-03 is that the lens is an M dwarf with mass 
of about 0.1-0.3 Mq located 1-5 kpc from the Earth and be- 
longing to the thick disk of the Galaxy, similar to the event 
MACHO-LMC-5, e.g. Drake et al. (2004), Gould (2004) and 
event EROS2-LMC-8 (Tisserand et al. 2007). 

With velocity of the order of 80 kms' 1 and location a 
few kpc from the Earth, the lens should have a noticeable 
proper motion of a few masyr -1 and therefore should be 
resolvable with the HST in a couple of years. This would 
give a good opportunity to confirm the nature of this event. 



4.3 OGLE-SMC-04 

This event was not previously detected neither by EWS (it 
occurred in the end of 2002 in the early season of the OGLE- 
III when the EWS was not yet operating), nor by any other 
surveys. Its light curve is relatively well covered by OGLE- 
III /-band data, but there are no observations taken during 
the event in V-band. 

The event was present in the EROS2 database, however 
only in their B-band. Therefore we performed a multi-band 
fit using OGLE /-band and EROS B-band data, which re- 
sulted in instrumental colour of the source of (B — I)s = 
0.82 ±0.04. This was then transformed to standard OGLE's 
V-I following Tisserand et al. (2007) and yielded (V-I)s = 
1.02 ±0.08. The error-bar includes the systematic error from 
the transformation. The position of the source on the CMD 
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Figure 8. Distributions of blending parameter (source flux frac- 
tion in baseline flux) for selected levels of stellar density derived 
for the simulated OGLE-III SMC images based on deep archival 
HST images. The distributions are shown for three magnitude 
bins. The density levels shown correspond to stellar densities 
of log(7V»/CCD chip) = (4.9,4.6,4.4,3.9) for very dense, dense, 
medium and sparse fields, respectively. 



Figure 9. Correction factor (CF) for the number of monitored 
stars as a function of the density of the stars on a single CCD 
chip of the OGLE-III template image. Shown are the CFs for 
three magnitude bins (B for bright, M for medium and F for 
faint) and for the entire magnitude range (A for all). The legend 
shows the expression for the curve fit to each of the data sets. 



(Fig. 2) indicates it belongs to the Red Clump giant popula- 
tion of the SMC. 

The time-scale of the event is relatively short (£e = 
18.3±1.8), which favours a self-lensing (SL) scenario. Unless 
the lens is nearby (Dl <20 kpc), a typical time-scale of halo 
lensing events with masses around 0.4 Mq and sources in 
the SMC are well above 20 days. In SL the time-scales are 
usually a few times shorter owing to much larger projected 
velocity, vsl ~ 2000 km/s, compared to Vhaio ~ 200 km/s 
for halo lenses (Boutreux & Gould 1996). This places this 
event as a good candidate for a self-lensing microlensing 
interpretation. 

Another possibility if that the lens is a thick disk red 
dwarf given the fact the blending object must be redder that 
the Red Clump of the SMC. High-resolution imaging could 
potentially reveal the actual lens as the event happened al- 
ready almost 10 years ago giving enough time for lens and 
source to separate. 



5 BLENDING AND DETECTION 
EFFICIENCY 

As in the previous parts of the study of the OGLE microlens- 
ing data towards the Magellanic Clouds we carefully tackled 
the issue of blending in these crowded fields. Here we applied 
the method developed for Paper III. In brief, we simulated 
OGLE-III images for a range of stellar densities using deep 
luminosity functions of three representative fields from the 
HST Local Group Stellar Photometry archive (Holtzman 
et al. 2006). Locations of the selected fields are marked on 
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Figure 10. Microlensing events detection efficiency for example 
dense (solid lines, field SMC106.3) and sparse (dashed lines, field 
SMC111.7) fields of the SMC as observed by OGLE-III. Shown 
are curves for All Stars Sample and Bright Stars Sample (thin 
and thick lines, respectively). 

Fig. 1. Then we compared simulated OGLE images with the 
input catalogue of stars and were able to assign a number 
of real HST stars to each blended object detected on a sim- 
ulated frame. For each density level and in three magnitude 
bins, 14 < / < 17.5 (bright), 17.5 < I < 19 (medium) and 
19 < I < 21 (faint) we derived a distribution of flux ratio 
between each HST component and the blended object it was 
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residing in. The blending distributions are shown in Fig. 8. 
We were also able to calculate a mean ratio between the 
number of real stars to all objects composed of these stars. 
This ratio was then used as a correction factor (CF) for ob- 
served OGLE template images in order to estimate the real 
number of monitored stars in each field (see Table 1). The 
CF is shown in Fig. 9 for all stars (A) and three magnitude 
bins: bright (A), medium (M) and faint (F). 

With these blending distributions at hand we were able 
to perform the detection efficiency determination. For se- 
lected fields representing different density levels and for 
time-scales in range between 1 and 1000 days we simu- 
lated numerous microlensing events with the parameters be- 
ing drawn from realistic distributions. For the lensed flux 
we used the measurements of randomly selected star from 
the simulated field, which during the event was split be- 
tween blend and source according to the blending param- 
eter. This procedure also preserved any variability of the 
underlying flux, allowing for events with variable baseline 
(see e.g., Wyrzykowski et al. 2006). 

The efficiency curves for representative dense and sparse 
fields and for the All Stars and Bright Stars samples are 
shown in Fig. 10. The detection efficiencies do not change 
much with the density of the field due to relatively low gra- 
dient in the stellar density in the OGLE-III SMC fields. The 
calculated efficiency is somewhat lower when compared with 
the one obtained for OGLE-III LMC data, mainly because 
here we used a different magnitude cut on star selection for 
the All Stars Sample (21.0 mag, compared with 20.4 mag 
for the LMC). 



6 OPTICAL DEPTH 

Having the three events and their time-scales as well as the 
detection efficiency of each of them we are able to calculate 
the average optical depth towards the SMC based on OGLE- 
III data. In order to evaluate Eq. 2 we also need the esti- 
mated number of total monitored stars, derived in the pre- 
vious section, and the total time of observation, T o b s =2870 
days. 

First, we calculated the optical depth for the All Stars 
Sample (N, = 5 971 776, down to 7=21 mag) using the 
detection efficiencies obtained in our simulations. Because 
our simulations utilised /-band data solely we used £e ob- 
tained in 5-parameter fit for all three events. This led to 
total tsmc-oiii = 1.17 ± 0.91 x 10~ 7 . The error on r was 
calculated following Han & Gould (1995) and reflects only 
uncertainty caused by the small number of events used. 

Then, similarly as in Papers I, II and III, we cor- 
rected the efficiency for the fact our simulations were only 
considering single-point events and that lensing by a bi- 
nary (and other exotic effects) was not included. Efficien- 
cies were scaled down by 10 per cent and the optical depth 
for efficiency corrected for binary events was tsmc-oiii = 
1.30 ± 1.01 x 10~ 7 . Details of the calculations for each 
event are shown in Table 7. 

For the Bright Stars Sample (iV* = 1 702 724, down to 
7=19.3 mag) the optical depth calculated for two events 
(#2 and #4) is tsmc-oiii = 1-52 ± 1.40 x 10~ 7 and 
tsmc-oiii = 1.69 ± 1.55 x 10~ 7 , when efficiencies were cor- 
rected for lack of sensitivity for binary events. Within the 



Table 7. The optical depth for the three events found in the 
OGLE-III SMC data calculated for the All Stars Sample. The 
columns show event name, its time-scale, detection efficiency and 
individual event's contribution to the total optical depth. 



event *e € (*e) t« X 10 7 

efficiency not corrected for binary events 

OGLE-SMC-02t 195.6±1.9 0.234942 0.76 

OGLE-SMC-03t 45.5±6.2 0.155922 0.27 

OGLE-SMC-04 18.60t^ 0.123365 0.14 
total tsmc-oiii 1.17 ±0.91 

efficiency corrected for binary events 
OGLE-SMC-02t 195.6±1.9 0.211448 0.85 
OGLE-SMC-03t 45.5±6.2 0.140330 0.30 
OGLE-SMC-04 18.60ti;g5 0.111028 0.15 



total tsmc-oiii 1-30 ± 1.01 



1 Time-scale and efficiency for these events are taken from a stan- 
dard fit. 

(large) error-bars, the estimates of the optical depth for the 
All Stars and Bright Stars Samples are in agreement. How- 
ever, it is clearly seen that higher detection efficiency does 
not completely counterbalance the smaller number of stars 
and smaller number of events in the Bright Stars Sample. 

Note that the binary-lens event OGLE-SMC-02 was in- 
cluded in the optical depth calculations above as if it were a 
standard event. This approximation is justified because even 
though the light curve of this event deviates a little from the 
standard "Paczyhski" curve, the event passes through our 
automated pipeline, which does not consider any exotic ef- 
fects. Moreover, the time-scale obtained in the standard fit 
is very close to the one obtained in the full model by Dong 
et al. (2007). Therefore it was also possible to derive the 
detection efficiency for this event. The same applies to the 
parallax event OGLE-SMC-03 



7 DISCUSSION 

In the 8 years of OGLE-III data covering the Small Magel- 
lanic Cloud and its surroundings we detected three convinc- 
ing microlensing events candidates. 

This is definitively the best set of candidates presented 
in the series of our papers concerning microlensing searches 
towards the Magellanic Clouds with OGLE data. The mi- 
crolensing nature of all three events is very difficult to dis- 
prove - none of their sources is located in the region of poten- 
tial contaminants "Blue Bumpers" (Fig. 2) and their light 
curves covering at least 8 years (12 years for OGLE-SMC-02 
as it was also monitored during the OGLE-II) show no other 
additional bumps. Additional data spanning for a couple of 
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more years available for all three events from the MACHO 
and EROS groups also show no further bumps. 

However, still the most difficult question to answer is 
where these events originated from. In principle, because 
standard microlensing model can not tell us where the source 
and the lens are located, there are a number of different 
combinations possible, including the SMC, the halo of the 
Milky Way and the disk of the Milky Way. In the history 
of the Magellanic Clouds microlensing there were already 
examples of confirmed cases of each of the combinations, 
e.g., source and lens from the SMC (Assef et al. 2006), disk 
lens towards the LMC (Kallivayalil et al. 2004). 

In the case of our events, the new one, OGLE-SMC-04, 
is probably the best candidate for self-lensing event given 
its time-scale and location in the Red Clump stars of the 
SMC, however, we can not exclude the thick disk red dwarf 
scenario or MACHO. The parallax effect detected in the 
event OGLE-SMC-03 help constrain the lens to most likely 
be a typical M-dwarf star from the Galactic disk, however a 
dark lens can not be entirely ruled out. Finally, even though 
OGLE-SMC-02 is the likely to be a halo lens, Dong et al. 
(2007), who modelled it in detail, did not exclude the self- 
lensing option completely. 

The effect of self-lensing for the SMC is not yet well un- 
derstood. However, due to the alignment of the SMC along 
the line-of-sight the contribution of SL to r is expected to be 
higher than in the face-on LMC and be at least 1.0xl0~ 7 , 
in average (Palanque-Delabrouille et al. 1998, Graff & Gar- 
diner 1999). This is in good agreement with the value of 
the optical depth we measured for all three events found in 
the OGLE-III data. Therefore, we expect the signal from 
the self-lensing or, more generally, the background of lens- 
ing by known stellar populations, is close to what we have 
detected. However, because the nature of our events is not 
yet firmly confirmed we can derive an upper limit on dark 
matter compact halo objects, following Alcock et al. (2001), 
Tisserand et al. (2007) and Papers I, II and III. 

Assuming the MACHOs mass distribution function 
from model "S" from Alcock et al. (2000) and using the 
mean detection efficiency for our SMC fields we derived the 
number of expected events due to MACHOs as a function of 
their mass. The mass in such events can be approximately 
translated to a time-scale as log M = 2 log(te/70). The num- 
ber of expected events compared with the observed (or not 
observed) signal can be converted to an upper limit on the 
MACHO 's halo mass fraction using enhanced Poisson statis- 
tics with background signal (Feldman & Cousins 1998). The 
derived upper limit for OGLE-III SMC data is shown in Fig. 
11 and was calculated assuming we expect about 3 events 
from the self-lensing/background. For a "typical" MACHO 
with mass of 0.4 Mq the fraction is less than 37 per cent, 
but for lower masses, between 0.01 and 0.1 Mq it is below 20 
per cent. Figure 11 shows also an upper limit on MACHOs 
derived by EROS (Tisserand et al. 2007) and MACHO sig- 
nal claimed by the MACHO group at 95 per cent confidence. 
For the latter we show both the original value from Bennett 
(2005) (for t — 1.0 i 0.3 x 10~ 7 ) for all MACHO events be- 
ing due to MACHOs and the value corrected following the 
rejection of one of the original MACHO events (MACHO- 
LMC-7, see Paper III). For the corrected value, the MACHO 
fraction in the halo becomes 18 per cent, compared to the 
original 20 per cent. If the background lensing (Galactic disk 



and LMC self-lensing) was extracted from the MACHO sig- 
nal (at a conservative level suggested in Bennett 2005 of 
0.24 x 10~ 7 ) and the compromised event was rejected, the 
MACHO group data can be interpreted that only about 14 
per cent of MACHOs may compose the Galactic halo. 

We can obtain a much tighter limit from OGLE by com- 
bining all OGLE-II and OGLE-III results for both Large and 
Small Magellanic Clouds. For each set of data we computed 
the expected number of events for a given mass of the deflec- 
tor and combined it with about 8 events expected to be due 
to the background lensing populations. It yielded / < 0.06 
(at 95 per cent C.L.) for masses around 0.4 Mq and reached 
the minimum of / < 0.04 at mass range between 0.01 and 
0.15 Mq. Figure 12 shows the combined OGLE limit along 
with the corrected MACHO group signal and EROS upper 
limit with the both axes zoomed in at the intersection of 
results from all surveys. This result, therefore, pushes the 
upper limit on MACHO mass fraction in the halo down to 
a value similar to the one derived by EROS. The OGLE re- 
sult is a somewhat more sensitive to MACHOs in the higher 
mass regime (between 1 and 10 Mq), and is significantly less 
sensitive for masses below 0.01 Mq. 

When looking at all events in more detail, in the most 
likely case, we can safely attribute to non-MACHO lens- 
ing (self-lensing or disk lensing) all the events detected in 
the LMC and SMC during OGLE-II and OGLE-III except 
OGLE-SMC-02. In such case, the optical depth due to this 
single dark matter event would be t$mc#2 ~ 0.12 x 10~ 7 . 
If the MW halo were composed only from such lenses with 
mass of around 10 Mq, such an optical depth would mean 
they contribute no more than 2 per cent of the total mass 
of the dark halo. This is in agreement with the expected 
mass fraction of black-holes in the total mass budget of the 
Galaxy varying from around 5 per cent (Sartore & Treves 
2010) for all massive stellar remnants to 0.4 per cent for a 
standard initial mass function (e.g., Bastian et al. 2010) in- 
tegrated above 6 Mq. It is also a much tighter limit than 
the limit of 40 per cent derived by the MACHO group on 
objects below 30 Mq (Alcock et al. 2001), obtained based on 
a lack of long-term events. The value obtained for a single 
microlensing event, however, is only crudely estimated, and 
should be assessed more elaborately, e.g., with simulations of 
long events in the OGLE-II and OGLE-III data combined. 
In practice, more events with long time-scales are needed in 
order to put better and more statistically sound constraints 
on the massive stellar remnants abundance in the halo. Con- 
tinued observations of the Magellanic Clouds by MOA and 
OGLE-IV surveys will hopefully yield more potential black- 
hole events residing in the halo if they are common enough. 



8 CONCLUSIONS 

We analysed OGLE-III data covering the Small Magellanic 
Cloud and found three convincing candidates for microlens- 
ing events, among which one was not known before (OGLE- 
SMC-04). Two of the events are likely to be caused by self- 
lensing or have lenses located in the Galactic disk. OGLE- 
SMC-03 is actually the first event towards the Small Mag- 
ellanic Cloud exhibiting clear parallax signal, which makes 
the lens likely to be in the thick disk. The third one (OGLE- 
SMC-02) is a potential candidate for a binary black hole lens 
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Figure 11. Fraction of the mass of MACHOs in the halo as derived from OGLE-III SMC data. Black solid curves show an upper limit 
from the OGLE-III SMC data assuming the background self-lensing signal of three events. Also shown is the upper limit from EROS 
(blue dashed curve), original measurement by MACHO (small cross and dotted curve) and the MACHO result corrected for the fact 
that one of their events was rejected by OGLE-III data (dashed curve, big cross). 
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Figure 12. Inclusion region for the fraction of the mass of MACHOs in the halo for combined OGLE-II and OGLE-III data for LMC 
and SMC (solid curve, pink/dark blue) when all OGLE events are attributed to expected self-lensing/background signal. The EROS 
upper limit and the MACHO signal are shown as in Fig. 11 in light/dark blue and green, respectively. Also shown is the fraction of mass 
due to BH-candidate event OGLE-SMC-02. 
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from the halo, but its self-lensing origin can not be totally 
excluded. 

Following a detailed analysis of the blending in the SMC 
we derived the optical depth for the All Stars Sample of 
tsmc-oiii = 1-30 ± 1.01 x 1CP 7 , which is in agreement with 
expected self-lensing signal towards the SMC. The upper 
limit on the dark matter in form of MACHOs derived for 
the OGLE-III SMC data alone was about 20 per cent for 
masses below 0.1 Mq. 

In this paper we concluded the studies of the microlens- 
ing searches conducted so far by the OGLE project towards 
the Magellanic Clouds. We presented the final result com- 
bining all available data and deriving new constraints on the 
fraction of massive compact halo objects in the Galactic halo 
of 6 per cent for M = 0.1 — 0.4 Mq and below 4 per cent for 
lower masses. For MACHOs with mass of 1 M and 20 M© 
the upper limits are / < 9 and / < 20 per cent, respectively. 
13 years of OGLE observations indicate the MACHO halo 
fraction is well below the value suggested by MACHO group 
(14 per cent after correction due to rejection of one of their 
events owing to a second bump for the same events found 
in the OGLE data) and in agreement with the limit derived 
by EROS survey. 

Our result indicates that baryonic dark matter in the 
form of relics of stars and very faint objects in the sub-solar 
mass regime is unlikely to inhabit the Milky Way's dark 
matter halo in any significant numbers. The presence of the 
black hole lens candidate towards the SMC agrees with the 
expected < 2 per cent contribution of black holes to the 
mass of the Galactic halo. 

With the OGLE project now in its fourth phase we 
hope the sensitivity to extremely long events will improve 
significantly within the next years when combined with the 
historic OGLE data. Long duration events will also be easily 
detectable in near-real-time by Gaia satellite, due to launch 
in 2013. Detailed follow-up observations of the alerts trig- 
gered by OGLE and Gaia will be crucial for revealing the 
true nature of detected events. This should result in an in- 
crease in the number of potential black-hole lenses or allow 
us to rule out heavy dark matter compact objects as well 
and close that topic definitively. 
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